Characterization of Higher-order Chromatin Structure in Bone Differentiation and Breast Cancer: A Dissertation

Higher-order genome organization is important for the regulation of gene expression by bringing different cis-regulatory elements and promoters in proximity. The establishment and maintenance of long-range chromatin interactions occur in response to cellular and environmental cues with the binding of transcription factors and chromatin modifiers. Understanding the organization of the nucleus in differentiation and cancer has been a long standing challenge and is still not well-understood. In this thesis, I explore the dynamic changes in the higher-order chromatin structure in bone differentiation and breast cancer. First, we show dynamic chromatin contact between a distal regulatory element and the promoter of Runx2 gene, which encodes the Runtrelated transcription factor 2 (RUNX2) that is essential for bone development. Next, via using a genome-wide approach, we show that breast cancer cells have altered long-range chromatin contacts among small, gene-rich chromosomes and at telomeres when compared with mammary epithelial cells. Furthermore, we assess the changes in nuclear structure and gene expression of breast cancer cells following Runt-related transcription factor 1 (RUNX1) deficiency, an event frequently observed in breast cancer. Finally, I present the role of the central ATPase subunit of the SWI/SNF complex, SMARCA4 (BRG1), in mediating nuclear structure and gene expression. Taken together, the research presented in this thesis reveals novel insight and paradigm for the dynamic changes in disease and differentiation, as well as uncovers previously unidentified roles for two chromatin regulatory proteins, RUNX1 and SMARCA4

The PPARgamma locus makes long-range chromatin interactions with selected tissue-specific gene loci during adipocyte differentiation in a protein kinase A dependent manner

Differentiation signaling results in reprogramming of cellular gene expression that leads to morphological changes and functional specialization of a precursor cell. This global change in gene expression involves temporal regulation of differentiation-specific genes that are located throughout the genome, raising the idea that genome structure may also be re-organized during cell differentiation to facilitate regulated gene expression. Using in vitro adipocyte differentiation as a model, we explored whether gene organization within the nucleus is altered upon exposure of precursor cells to signaling molecules that induce adipogenesis. The peroxisome proliferator-activated receptor gamma (PPARgamma) nuclear hormone receptor is a master determinant of adipogenesis and is required for adipose differentiation. We utilized the chromosome conformation capture (3C) assay to determine whether the position of the PPARgamma locus relative to other adipogenic genes is changed during differentiation. We report that the PPARgamma2 promoter is transiently positioned in proximity to the promoters of genes encoding adipokines and lipid droplet associated proteins at 6 hours post-differentiation, a time that precedes expression of any of these genes. In contrast, the PPARgamma2 promoter was not in proximity to the EF1alpha promoter, which drives expression of a constitutively active, housekeeping gene that encodes a translation elongation factor, nor was the PPARgamma2 promoter in proximity to the promoter driving the expression of the C/EBPalpha regulatory protein. The formation of the long-range, intergenic interactions involving the PPARgamma2 promoter required the regulatory factor C/EBPbeta, elevated cyclic AMP (cAMP) levels, and protein kinase A (PKA) signaling. We conclude that genome organization is dynamically remodeled in response to adipogenic signaling, and we speculate that these transient inter-genic interactions may be formed for the purposes of selecting some of the transcriptionally silent tissue-specific loci for subsequent transcriptional activation

The connection between BRG1, CTCF and topoisomerases at TAD boundaries

The eukaryotic genome is partitioned into topologically associating domains (TADs). Despite recent advances characterizing TADs and TAD boundaries, the organization of these structures is an important dimension of genome architecture and function that is not well understood. Recently, we demonstrated that knockdown of BRG1, an ATPase driving the chromatin remodeling activity of mammalian SWI/SNF enzymes, globally alters long-range genomic interactions and results in a reduction of TAD boundary strength. We provided evidence suggesting that this effect may be due to BRG1 affecting nucleosome occupancy around CTCF sites present at TAD boundaries. In this review, we elaborate on our findings and speculate that BRG1 may contribute to the regulation of the structural and functional properties of chromatin at TAD boundaries by affecting the function or the recruitment of CTCF and DNA topoisomerase complexes

A TAD boundary is preserved upon deletion of the CTCF-rich Firre locus

The binding of the transcriptional regulator CTCF to the genome has been implicated in the formation of topologically associated domains (TADs). However, the general mechanisms of folding the genome into TADs are not fully understood. Here we test the effects of deleting a CTCF-rich locus on TAD boundary formation. Using genome-wide chromosome conformation capture (Hi-C), we focus on one TAD boundary on chromosome X harboring ~ 15 CTCF binding sites and located at the long non-coding RNA (lncRNA) locus Firre. Specifically, this TAD boundary is invariant across evolution, tissues, and temporal dynamics of X-chromosome inactivation. We demonstrate that neither the deletion of this locus nor the ectopic insertion of Firre cDNA or its ectopic expression are sufficient to alter TADs in a sex-specific or allele-specific manner. In contrast, Firre’s deletion disrupts the chromatin super-loop formation of the inactive X-chromosome. Collectively, our findings suggest that apart from CTCF binding, additional mechanisms may play roles in establishing TAD boundary formation

Cancer-testis gene expression is associated with the methylenetetrahydrofolate reductase 677 C\u3eT polymorphism in non-small cell lung carcinoma

BACKGROUND: Tumor-specific, coordinate expression of cancer-testis (CT) genes, mapping to the X chromosome, is observed in more than 60% of non-small cell lung cancer (NSCLC) patients. Although CT gene expression has been unequivocally related to DNA demethylation of promoter regions, the underlying mechanism leading to loss of promoter methylation remains elusive. Polymorphisms of enzymes within the 1-carbon pathway have been shown to affect S-adenosyl methionine (SAM) production, which is the sole methyl donor in the cell. Allelic variants of several enzymes within this pathway have been associated with altered SAM levels either directly, or indirectly as reflected by altered levels of SAH and Homocysteine levels, and altered levels of DNA methylation. We, therefore, asked whether the five most commonly occurring polymorphisms in four of the enzymes in the 1-carbon pathway associated with CT gene expression status in patients with NSCLC. METHODS: Fifty patients among a cohort of 763 with NSCLC were selected based on CT gene expression status and typed for five polymorphisms in four genes known to affect SAM generation by allele specific q-PCR and RFLP. RESULTS: We identified a significant association between CT gene expression and the MTHFR 677 CC genotype, as well as the C allele of the SNP, in this cohort of patients. Multivariate analysis revealed that the genotype and allele strongly associate with CT gene expression, independent of potential confounders. CONCLUSIONS: Although CT gene expression is associated with DNA demethylation, in NSCLC, our data suggests this is unlikely to be the result of decreased MTHFR function

Epigenetic control of cell cycle-dependent histone gene expression is a principal component of the abbreviated pluripotent cell cycle

Self-renewal of human pluripotent embryonic stem cells proceeds via an abbreviated cell cycle with a shortened G(1) phase. We examined which genes are modulated in this abbreviated period and the epigenetic mechanisms that control their expression. Accelerated upregulation of genes encoding histone proteins that support DNA replication is the most prominent gene regulatory program at the G(1)/S-phase transition in pluripotent cells. Expedited expression of histone genes is mediated by a unique chromatin architecture reflected by major nuclease hypersensitive sites, atypical distribution of epigenetic histone marks, and a region devoid of histone octamers. We observed remarkable differences in chromatin structure--hypersensitivity and histone protein modifications--between human embryonic stem (hES) and normal diploid cells. Cell cycle-dependent transcription factor binding permits dynamic three-dimensional interactions between transcript initiating and processing factors at 5\u27 and 3\u27 regions of the gene. Thus, progression through the abbreviated G(1) phase involves cell cycle stage-specific chromatin-remodeling events and rapid assembly of subnuclear microenvironments that activate histone gene transcription to promote nucleosomal packaging of newly replicated DNA during stem cell renewal

The bone-specific Runx2-P1 promoter displays conserved three-dimensional chromatin structure with the syntenic Supt3h promoter

Three-dimensional organization of chromatin is fundamental for transcriptional regulation. Tissue-specific transcriptional programs are orchestrated by transcription factors and epigenetic regulators. The RUNX2 transcription factor is required for differentiation of precursor cells into mature osteoblasts. Although organization and control of the bone-specific Runx2-P1 promoter have been studied extensively, long-range regulation has not been explored. In this study, we investigated higher-order organization of the Runx2-P1 promoter during osteoblast differentiation. Mining the ENCODE database revealed interactions between Runx2-P1 and Supt3h promoters in several non-mesenchymal human cell lines. Supt3h is a ubiquitously expressed gene located within the first intron of Runx2. These two genes show shared synteny across species from humans to sponges. Chromosome conformation capture analysis in the murine pre-osteoblastic MC3T3-E1 cell line revealed increased contact frequency between Runx2-P1 and Supt3h promoters during differentiation. This increase was accompanied by enhanced DNaseI hypersensitivity along with RUNX2 and CTCF binding at the Supt3h promoter. Furthermore, interplasmid-3C and luciferase reporter assays showed that the Supt3h promoter can modulate Runx2-P1 activity via direct association. Taken together, our data demonstrate physical proximity between Runx2-P1 and Supt3h promoters, consistent with their syntenic nature. Importantly, we identify the Supt3h promoter as a potential regulator of the bone-specific Runx2-P1 promoter. Acids Research

SMARCA4 regulates gene expression and higher-order chromatin structure in proliferating mammary epithelial cells

The packaging of DNA into chromatin plays an important role in transcriptional regulation and nuclear processes. Brahma-related gene-1 SMARCA4 (also known as BRG1), the essential ATPase subunit of the mammalian SWI/SNF chromatin remodeling complex, uses the energy from ATP hydrolysis to disrupt nucleosomes at target regions. Although the transcriptional role of SMARCA4 at gene promoters is well-studied, less is known about its role in higher-order genome organization. SMARCA4 knockdown in human mammary epithelial MCF-10A cells resulted in 176 up-regulated genes, including many related to lipid and calcium metabolism, and 1292 down-regulated genes, some of which encode extracellular matrix (ECM) components that can exert mechanical forces and affect nuclear structure. ChIP-seq analysis of SMARCA4 localization and SMARCA4-bound super-enhancers demonstrated extensive binding at intergenic regions. Furthermore, Hi-C analysis showed extensive SMARCA4-mediated alterations in higher-order genome organization at multiple resolutions. First, SMARCA4 knockdown resulted in clustering of intra- and inter-subtelomeric regions, demonstrating a novel role for SMARCA4 in telomere organization. SMARCA4 binding was enriched at topologically associating domain (TAD) boundaries, and SMARCA4 knockdown resulted in weakening of TAD boundary strength. Taken together, these findings provide a dynamic view of SMARCA4-dependent changes in higher-order chromatin organization and gene expression, identifying SMARCA4 as a novel component of chromatin organization

Chromatin interaction analysis reveals changes in small chromosome and telomere clustering between epithelial and breast cancer cells

BACKGROUND: Higher-order chromatin structure is often perturbed in cancer and other pathological states. Although several genetic and epigenetic differences have been charted between normal and breast cancer tissues, changes in higher-order chromatin organization during tumorigenesis have not been fully explored. To probe the differences in higher-order chromatin structure between mammary epithelial and breast cancer cells, we performed Hi-C analysis on MCF-10A mammary epithelial and MCF-7 breast cancer cell lines. RESULTS: Our studies reveal that the small, gene-rich chromosomes chr16 through chr22 in the MCF-7 breast cancer genome display decreased interaction frequency with each other compared to the inter-chromosomal interaction frequency in the MCF-10A epithelial cells. Interestingly, this finding is associated with a higher occurrence of open compartments on chr16-22 in MCF-7 cells. Pathway analysis of the MCF-7 up-regulated genes located in altered compartment regions on chr16-22 reveals pathways related to repression of WNT signaling. There are also differences in intra-chromosomal interactions between the cell lines; telomeric and sub-telomeric regions in the MCF-10A cells display more frequent interactions than are observed in the MCF-7 cells. CONCLUSIONS: We show evidence of an intricate relationship between chromosomal organization and gene expression between epithelial and breast cancer cells. Importantly, this work provides a genome-wide view of higher-order chromatin dynamics and a resource for studying higher-order chromatin interactions in two cell lines commonly used to study the progression of breast cancer

Enhancers in the Peril lincRNA locus regulate distant but not local genes

Abstract Background Recently, it has become clear that some promoters function as long-range regulators of gene expression. However, direct and quantitative assessment of enhancer activity at long intergenic noncoding RNA (lincRNA) or mRNA gene bodies has not been performed. To unbiasedly assess the enhancer capacity across lincRNA and mRNA loci, we performed a massively parallel reporter assay (MPRA) on six lincRNA loci and their closest protein-coding neighbors. Results For both gene classes, we find significantly more MPRA activity in promoter regions than in gene bodies. However, three lincRNA loci, Lincp21, LincEnc1, and Peril, and one mRNA locus, Morc2a, display significant enhancer activity within their gene bodies. We hypothesize that such peaks may mark long-range enhancers, and test this in vivo using RNA sequencing from a knockout mouse model and high-throughput chromosome conformation capture (Hi-C). We find that ablation of a high-activity MPRA peak in the Peril gene body leads to consistent dysregulation of Mccc1 and Exosc9 in the neighboring topologically associated domain (TAD). This occurs irrespective of Peril lincRNA expression, demonstrating this regulation is DNA-dependent. Hi-C confirms long-range contacts with the neighboring TAD, and these interactions are altered upon Peril knockout. Surprisingly, we do not observe consistent regulation of genes within the local TAD. Together, these data suggest a long-range enhancer-like function for the Peril gene body. Conclusions A multi-faceted approach combining high-throughput enhancer discovery with genetic models can connect enhancers to their gene targets and provides evidence of inter-TAD gene regulation